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Abstract 

We investigate the dark matter constraints for the nonminimal SUSY standard model 
(NMSSM). The cosmologically restricted mass spectra of the NMSSM are compared to 
the minimal SUSY standard model (MSSM). The differences of the two models concerning 
the neutralino, sfermion and Higgs sector are discussed. The dark matter condition leads 
to cosmologically allowed mass ranges for the SUSY particles in the NMSSM: m^o < 
300 GeV,mg^ < 300 Ge^, 300 GeV < m^^ < 1900 GeV, 200^61^ < m^-^ < 1500 GeV", 
350 Gel^ < rrig < 2100 Gel^ and for the mass of the lightest scalar Higgs ms^ < lAOGeV . 



i) Introduction. Supersymmetry (SUSY) is suggested to solve the hierarchy problem 
of the standard model (SM), if it is embedded in a grand unified theory (GUT). The 
simplest SUSY extension of the SM is the minimal SUSY standard model (MSSM) [0, 0. 
In the nonminimal SUSY standard model (NMSSM) the Higgs sector of the MSSM 
is extended by a Higgs singlet N. By introducing a Higgs singlet N the /i-parameter of 
the MSSM (/i-problem) can be dynamically generated via fi = Xx with the vaccum ex- 
pectation value (VEV) (N) = x and A being the Yukawa coupling of the Higgs fields. 
The phenomenology of the NMSSM has been analysed in several papers [Q - SUSY 
models with Higgs singlets can be derived from superstring inspired or SU{5) x f/(l) 
GUT models and they offer the possibility of spontaneous breaking of the CP symme- 
try. The discrete Z3 symmetry of the NMSSM causes a domain wall problem, because 
the Z3 symmetry is spontaneously broken during the electoweak phase transition in the 
early universe. The domain wall problem of the NMSSM and possible solutions to it are 



discussed in [14 



In addition to the particle spectrum of the MSSM there are two extra Higgses and one 
extra neutralino in the NMSSM, which can mix with the other Higgses and neutralinos, 
and thus modify their properties. In many SUSY models the lightest SUSY particle (LSP) 
is a neutralino, which is a good candidate for cold dark matter. In the NMSSM the LSP 
neutralino can have a larger singlino portion. The cosmology of these LSP singlinos has 
been discussed in [IT3 at a given fixed low energy scale. The experimentally and cosmo- 
logically allowed parameter space of the NMSSM and the dark matter neutralinos have 
been investigated in using RG evolutions. 

In the present paper we show, how the dark matter condition restricts the mass spec- 
tra of the SUSY particles in the NMSSM. The cosmologically restricted mass spectra of 
the MSSM and NMSSM are compared and differences are discussed. In contrast to the 
previous investigation on NMSSM dark matter |]16|, we search for solutions which are 
the (global) minimum of the effective 1 loop Higgs potential and whose parameters are 
evolved with the SUSY RGEs from the GUT scale to the electroweak scale using the lat- 
est experimental constraints. Furthermore we consider all possible decay channels which 
occur in neutralino annihilation. 
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ii) The NMSSM. The NMSSM § is a super symmetric SU{?,)c x SU{2)i x f/(l)y gauge 
theory with two Higgs doublets ifi, if 2 and a Higgs singlet A^. This model is defined by 
the superpotential (with only dimensionless couplings) 

W = hdH^ eQ i) -KH^ eQU + XH^ eH2N -\kN^ (1) 

(e is the antisymmetric tensor with ei2 = 1) and the SUSY soft breaking terms 

Csoft = -mj,jHi\''-mljH2\^-m%\N\^-M^\Q\^-Mlj\U\''-Ml\D\^ 
-hdAdHf eQD + K HjeQU 
+ XAxH^ eH2N + \kAkN^ 

+ i M3 A3 A3 + ^ M2 + ^ Ml Ai Ai + h.c. (2) 

The effective 1 loop Higgs potential consists of the tree-level potential and the radiative 
corrections 

Vlloop = Vtree + Vrad- (3) 

The contributions of the top quark and stops ti^2 are considered in the radiative corrections 
to the effective Higgs potential 

64 7r2 ^ ^ ' '"'Q 

i 

where the sum is taken over all particles and antiparticles with field-dependent mass rrii, 
spin Si and color degrees of freedom Cj. The electroweak gauge-symmetry SU{2)j x U (l)y 
is spontaneously broken to the electromagnetic gauge-symmetry U{l)em by the Higgs 
VEVs (H^) = Vi with i = 1,2 and (A^) = x. The three minimum conditions for the VEVs 
have the form: 

- m i' = 7^—^ X X 5 

2 ^ tan2/3-l ^ ' 

. ^ 2Xx{Ax + kx) 

"'""^^^ + E^ + mjj^ + E^ + 2 X^ x^ + X^ ^ ' 

{m% + S^) x^ - kAkX^ + 2k'^ x^ + X"^ x"^ v'^ -^{Ax + 2kx) Xx f ^ sin(2/3) = (7) 



Vrad = 771^ E ("l)'"" (25, + 1) In(^), (4) 
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where v = \/vf + V2 = 174 GeV, tan/3 = V2/V1 and = dVrad/dvf. 

The mass of the neutrahnos follows from the following part of the lagrangian 

1 



-^^ + h.c. 



(9) 



The symmetric mass matrix M of the neutralinos in the basis given in (9) has the form: 

/ Ml —mz sin 6w cos (3 mz sin 6w sin /5 \ 

M2 mz cos 6w cos /3 — cos 6w sin /5 

— sin 9w cos /5 cos 6w cos /5 \x A f 2 

sin 6^14/ sin /? — cos 6*1^ sin /? A a? A f 1 

\ \V2 \vi -2k X J 

With fi = X X the first 4x4 submatrix recovers the mass matrix of the MSSM p[. The 
mass of the neutralinos is here obtained by diagonalising the mass matrix M with the 
orthogonal matrix A^. (Then some mass eigenvalues may be negative |T^.) 



C 



1 ~o ~o 

2 Xi Xi 



(10) 



with x° = Nij^^j and Mdiag = N M . 



Gauge coupling unification in non-supersymmetric GUTs is already excluded by the pre- 
cision measurements at LEP |jl8|. In supersymmetric GUTs gauge coupling unification 
g^{Mx) = 95 ^ 0.72 is possible at a scale Mx ^ 1.6 x 10^^ GeV. For the SUSY soft 
breaking parameters we suppose universality at the GUT scale Mx'. 



mi{Mx) 
Ma{Mx] 
MMx) 



: mo 
= mi/2 
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{Ax{Mx) = -Ao). 

The Yukawa couplings at the GUT scale take the values A(Mx) = Aq, k{Mx) = and 
ht{Mx) = hto. With the SUSY RGEs ^ the SUSY soft breaking parameters and the 
Yukawa couplings are evolved from the GUT scale to the electroweak scale M^eak — 
100 GeV. At Mweak "we minimize the effective 1 loop Higgs potential. The 8 parameters 
of the NMSSM are mo, mi/2, Aq, Xq, ko, hto, tan/5 and x, with only 5 of them being 
independent, which in our procedure are taken as Ao, fco, hfQ, tan (3 and x. The remaining 3 



parameters (mo, mi/2, Aq) are numerically calculated from the three minimum conditions 
(5-7) following the method described in P: 

Since we take here tan/3 and x as (randomly chosen) fixed input parameters, the three 
minimum conditions (5-7) depend only on the low energy parameters rrii, Ai, A, k and ht 
generalized as pk{My,eak)- 



d Viioop 



FM = ^ (z = l-3). (13) 



dvi 

The low energy parameters A, k and ht are fixed by the SUSY RGEs and the remaining 
(randomly chosen) fixed input parameters Aq, k^ and hto- With the SUSY RGEs the low 
energy parameters nii and Ai are then only depending on the GUT scale parameters mo, 
mi/2 and Aq. Therefore all low energy parameters Pki^weak) can be written as functions 
fk (which can be constant) of the GUT scale parameters mo, mi/2 and Aq. 

Pk{Myjeak) = fk{mo,mii2,Ao). (14) 



The functions fk are determined by numerically solving the SUSY RGEs |jT9[. In the 
three minimum conditions (13) we then substitute the low energy parameters Pk{Mu,eak) 
by the functions fk and obtain the equations 

^i(/fe("^o,"^i/2,^o)) = G'i(mo,mi/2, Ao) = (i = 1 - 3). (15) 

The SUSY breaking parameters mo, mi/2 and Aq are then determined by numerically 
solving the equations Gj(mo, mi/2, ^o) =0 (i = 1 — 3). 

Because we assume ht ^ hb,hr, we restrict |tan/3 | to be smaller than 20 [^. For the 
singlet VEV x we choose the range |x| < 80 TeV 0, [11|. 

The physical minimum of the Higgs potential should be the global minimum. We check, 
whether the physical minimum is lower than the unphysical minima of Viioop(fi, 172, x) 
with at least one vanishing VEV M. 



The conditions |jT9[, that the minimum of the scalar potential does not break the conser- 
vation of charge and colour, can be formulated as follows: 

Al^ < 3(m^^ + m^^ + mfjj at scale Q ~ A^Jh^^ (16) 
A% < 3(m^j + mg^ + m^J at scale Q ~ Adjhd, (17) 
Al < 3{mjj^ + m\ + m|J at scale Q ~ AJK^. (18) 
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For the top quark pole mass we take the range 169 GeV < rrit < 181 GeV corresponding 
to the recent world average measurement The SUSY particles in the NMSSM have 
to fulfill the following experimental conditions: mj> > 41.8 GeV ^ij, rrig > 65 GeV p^ , 



rrig > 176 GeV (if m-g < 300 GeV) ^ and m,- > 70 GeV (for all 



mi > 58.3 GeV 0, > 173 GeV > 71.3 GeV mH+ > 43.5 GeV [0, 

fnsi > 70.7 GeV or the coupling of the lightest scalar Higgs Si to the Z boson is 
reduced compared to the SM ^ assuming visible decays with branching ratios like 
the SM Higgs 0, and 

^r(Z-.x°x;)<30Me\/ 



r(Z^X?X?)<7Mey |2 
BR{Z < 10-' , ^ (1, 1) m 



iii) i^arA; matter constraints. Infiationary cosmological models suggest a fiat universe 
with Q = p/pcrit = 1- Big bang nucleosynthesis restricts the baryon density to fltm-yons < 
0.1. The missing matter is called dark matter. One of the favoured theories to explain the 
structure formation of the universe is the cold + hot dark matter model (CHDM) |29| . 
In this model the dark matter consists of hot dark matter (massiv neutrinos) and cold 
dark matter (neutralinos) with Q^ot — 0.3 and ^Icoid — 0.65. With the Hubble constant 
Hq = 100 ho km s~^ Mpc~^ and the range 0.4 < ho < 1, the condition for the cosmic 
density of the neutralinos in the CHDM model reads 



0.1 < n^hl < 0.65. 



(19) 



With the annihilation cross section of the neutralinos (JannV = a + bv'^ {v is the relative 
velocity of the neutralinos) the cosmic neutralino density |]30| , |3ll , |32| can be calculated 

^^Mp{a + ^^)GeV' ^ ' 

where Mp = 1.221 x 10^^ GeV is the Planck mass and g^, ~ 81 is the effective number of 
degrees of freedom at Tfr- The freeze-out temperature of the neutralinos Tjr follows from 



the equation 30, 31, 32 



Xfr 



In I 0.0764 \m^o\ Mp y c (2 + c) 



•\/ 9* fr 



(21) 



with Xfr = \m^o\/Tfr and c = 1/2. 

In the MSSM the formulas for the annihilation cross section of the neutralinos into the 



different decay products X and Y can be found in |^ 



Crann(X^ ^ X Y)v = - 



\ACSo)\' + li\ACPo)\' + \ACP^)\' + \ACP2)\') 



(22) 



s ~ 4m? is the center of mass energy squared. S is a symmetry factor, which is 2 if 
X = Y. The partial wave amplitude A describes annihilation from an initial state "^^^^Lj 
with spin S, orbital angular momentum L and total angular momentum J. 
In the NMSSM the annihilation cross section of the neutralinos changes because of more 
particles in this model, which can mix, and modified vertices. In our numerical analysis we 
consider all possible decay channels. The relevant formulas for the NMSSM are obtained 
by substituting the MSSM couplings in |^ by the corresponding NMSSM couplings and 
by slightly modifying the partial wave amplitudes in |3^ and will be given in [Q. Here 



we only give the NMSSM couplings for the dominant decay channel 

Xi Xi ~^ fa fa- 

The fermions are produced by Z boson, scalar Higgs 5*4 (i = 1 — 3) and pseudoscalar 
Higgs Pa (a = 1, 2) exchange in the s-channel and sfermion /i 2 exchange in the t- and 
u-channel. When the LSP neutralino, which is denoted with index in instead of 1, is 



predominantly a bino (A^oi^ ~ 1 in eq. (11)), the annihilation cross section is dominated 
by sfermion exchange. The couplings of the unmixed sfermion /j to the fermion and 
LSP neutralino, are given by 

X„o = -V2g2 [T^a N02 - tan 9w (Tg, - e^J iVoi] (24) 
Yao = V2g2 tan 9w e /„ A^oi (25) 



.„o = 7^^^ ^04 (26) 



\/2 sin P rriw 
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Z^o = - ''"^^ iVo3 (27) 
V / cos p mw 

and refer to eqs. (A32a)-(A32c) in These couplings have the same form in the MSSM 
and NMSSM. The only difference follows from the matrix A^, which diagonalizes the mass 
matrix of the neutralinos. The couplings X^g, W'^^, Z'^q and Y^q of the mixed sfermion 
/i 2 to the fermion fa and LSP neutralino are easily given in terms of the unmixed case as 
described by eq. (A31) in [^. When the Higgsino portion of the LSP neutralino is larger 
(larger A'qs^, Nq^^ in eq. (11)), Z boson exchange becomes important for the annihilation 
cross section. The coupling of the Z boson to the LSP neutrahnos O^'^ (eq. (A5e) in |0l) 
has the same form in the MSSM and NMSSM, but a different matrix A^ enters: 

0'^^ = -In,,' + ^No^. (28) 

In the NMSSM the Higgs-couplings to the fermions have to be modified compared to the 
MSSM, 

2 sm p rriw 2 cos p 

, 92 rriu UI2 92 mg 

h'S.u = ■ n '^'S^d = ^ ■ l3Uj 

2 sm p iTLw 2 cos p rriuj 

The 3x3 mass matrix of the scalar Higgses Si [Q] in the basis {H^n-, -f^2ij) ^r) "with 
H^^ = \/2 Re{Hi} (i = 1-3) is diagonalised by the matrix . The 3x3 mass matrix of 
the pseudoscalar Higgses Pa and the Goldstone boson in the basis {H^j, H2J, Nj) with 
Hfj = ^/2Im{Hf} (i = 1-3) is here diagonahsed by the 3 x 3 matrix . These NMSSM 
couplings have to be substituted for the MSSM couplings given in eqs. (A33a)-(A33c) of 
113011 . The coupling of the scalar Higgs Sa to the neutralinos is different in the NMSSM 
and the MSSM. In the NMSSM this coupling can be written as 

Ts.^J = -U!,Q'l^ + U^S:^ + U!^z:^ (31) 
Q'l^ = [Ar,3 {g^ N,2 - 9v N,i) + V2 X N,, N,, + {i ^ j)] (32) 

^ 92 

S':^ = [Ar,4 {g^ N,2 - 9y iV,i) -V2XN,s N,, + {^ ^ j)] (33) 

^ 92 
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Z1, = V2 A iV,3 iV,4 + V2kNi^ iV,5 + j)l (34) 

^92 



to be substituted for the MSSM couplings given in eqs. (A9b,A9c) of [^. The couphng of 
the pseudoscalar Higgs Pa to the neutrahnos also changes. In the NMSSM this coupling 
has the following form 

Tp^ij = -Uai Qi'j + Ua2 -U^s Z-j (35) 
= 77^ [N^3 {92 N^2 - 9y N,i) - V2 \ N,, Nj, + (36) 

^ 92 

Si- = ^[N^i {92 Nj2 - 9y N,i) + V2XN,s N,, + (z ^ j)] (37) 

^92 

= z;;, (38) 

which has to be used instead of the MSSM coupling in eq. (A9a) of pO| . 

iv) Particle spectra. To obtain the mass spectra of the SUSY particles in the NMSSM, 
we randomly generate ~ 5.5x10^ points in the 5 dimensional parameter space of the input 
parameters Aq, hto, tan (3 and x. The SUSY breaking parameters mg, 1711/2 and Aq are 
then determined by numerically solving the three minimum equations (15) as described 
there. Having found a solution to the SUSY RGEs, whose effective 1 loop Higgs potential 
has a local minimum for the input parameters tan/3 and x, we calculate the masses of 
all SUSY particles 0, Q and Higgses [|, ||, |ll| and impose the additional theoretical 



and experimental constraints described above in ii). With all these requirements, about 
4900 solutions remain, from which about 2000 solutions are cosmologically acceptable. 



For the constrained MSSM in |^ a scatter-plot of solutions in the {mi/2,mQ) plane is 
shown. Upper bounds of 1.1 TeV for mi/2 and 600 GeV for mo are found. For mi/2 ~ 
100 GeV solutions with much larger mo are allowed due to Z-pole enhanced neutralino 
pair annihilation (m^o ^ m^/2). The comparable scatter-plot for the NMSSM (more 
details in [^]) shows only a few cosmologically allowed solutions with mo ^ Tni/2. The 
reason, why the case mo ^ rni/2 is somewhat suppressed in the NMSSM, is discussed in 
1^. It follows from the condition, that the minimum of the scalar potential does not break 
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the conservation of charge and colour, and from the condition that the physical minimum 
is lower than the symmetric vacuum, Aq > Qttiq. 

As a few representative examples for the importance of the dark matter constraints we 
show in the following some scatter-plots for the mass of SUSY particles. In the upper left 
picture of fig.0 the mass of the LSP neutralino versus the mass of the lighter selectron is 
shown for the MSSM. In the MSSM the LSP neutralino is mostly a bino. In the upper right 
picture the dark matter condition is imposed and gives un upper bound of 280 GeV for the 
LSP neutralino and a cosmologically allowed mass range of 100 GeV < rrie^ < 300 GeV 
for the lighter selectron. These bounds for the LSP neutralino and lighter selectron can 
be understood from the fact [^, that for a bino-like LSP neutralino the annihilation 
cross section is dominated by sfermion exchange. Heavier selectrons are cosmologically 
allowed, if Z-pole or Higgs-pole enhanced neutralino pair annihilation |^5[ occurs. 



For the NMSSM the lower left picture of fig.|l] shows the mass of the LSP neutralino 
versus the mass of the lighter selectron. Compared to the MSSM there are not so many 
solutions in the NMSSM with small LSP neutralino mass and large selectron mass, which 
follows from the suppression of mo ^ Tni/2 in the NMSSM §]. In the NMSSM some LSP 
neutralinos can have a larger singlino portion. Most of these singlinos are decoupled pure 
singlinos |TT|, but a smaller fraction of the singlinos can be mixed states (|T^, more 
details in ||3^). The LSP neutralinos with mass below 40 GeV and mg^ ~ 70 GeV in the 
lower left picture of fig.0 are lighter singlinos. 

The lower right picture of fig.|l] shows the NMSSM solutions, which survive the dark 
matter condition. The cosmologically allowed solutions in the MSSM and NMSSM look 
quite similar except that heavier selectrons are allowed in the former. With the dark 
matter condition nearly all decoupled pure LSP singlinos can be excluded (more details 



in ||3^), what can be seen for the lighter singlinos {m^o < AO GeV with mg^ ~ 70 GeV) 
by comparing the lower pictures of fig.[I|. 

Another difference between the MSSM and NMSSM is the enhanced lower bound of 
50 GeV for most of the LSP neutralinos in the NMSSM. The dark matter condition in 
the NMSSM in fig.|I] exludes nearly all LSP neutralinos with mass m^o ^ rriz/^, even if 
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they are not singlinos. The reason for this difference is the larger couphng Oqq in eq. (28) 
of the Z boson to the LSP neutrahnos in the NMSSM corresponding to neutrahnos with 
larger Higgsino portion. In the MSSM the coupling Oqq can be very small for nearly pure 
bino-like LSP neutrahnos, so that the neutralino pair annihilation becomes not too large 
from Z exchange near the Z-pole. The difference follows from the circumstance, that the 
NMSSM can be obtained from the MSSM in the limit A, — > 0, | x | ^ oo with Ax and kx 
fixed. In this limit the first two minimum conditions (5), (6) of the NMSSM correspond 
to the two minimum conditions of the MSSM, whereas the third minimum condition (7) 
of the NMSSM would be an extra constraint for the MSSM |]^, |TT[|. The third minimum 
condition restricts the parameter region of the NMSSM compared to the MSSM in such a 
way, that the coupling Oqq of the Z boson to the LSP neutrahnos is larger in the NMSSM 
for m^o ^ mzl'i- 

In the left picture of fig.0 we show the mass of the lighter u-squark and the lighter stop 
versus the mass of the gluino for the NMSSM. In contrast to the MSSM [^, the squark 
masses in the NMSSM are nearly proportional to the gluino mass. The already mentioned 
reason is, that in the NMSSM solutions with mo ^ ^1/2 are somewhat disfavoured. In 
the right picture of fig.0 the cosmologically allowed lighter u-squark and lighter stop 
masses are shown. In contrast to the MSSM, the dark matter condition improves the 
lower bound for the gluino to 350 GeV in the NMSSM. The reason is the enhanced lower 
bound of 50 GeV for most of the LSP neutrahnos in fig.0. Furthermore, the Higgsino 
or singlino portion of the lightest cosmologically allowed LSP neutrahnos is larger in the 
NMSSM than in the MSSM as already discussed. For these LSP neutrahnos with larger 
Higgsino or singlino portion, the gaugino mass M\ is larger than their mass m^o . A larger 
gaugino mass Mi leads to a larger gluino mass M3 in the NMSSM. Because the mass of 
the squarks are nearly proportional to the gluino mass in the NMSSM, the lower bound 
for the lighter u-squark improves to 300 GeV and the lower bound for the lighter stop 
improves to 200 GeV . This is in contrast to the MSSM, where the dark matter condition 
does not change the lower bounds. The cosmologically upper bound of 1900 GeV for the 
lighter u-squark, 1500GeV^ for the lighter stop and 2100GeV^ for the gluino in the NMSSM 
are similar to the corresponding upper bounds in the MSSM and are connected with the 



10 



cosmologically upper bound of 300 GeV for the lighter selectron and the bino-hke LSP 
neutrahno. 



In the left picture of fig.^ the mass of the lightest scalar Higgs 5*1 versus the mass of the 
gluino is shown for the NMSSM. The lightest scalar Higgs can be lighter than the MSSM 



Higgs bound of 62.5 GeV |^. In this case the lightest scalar Higgs is predominantly a 
Higgs singlet ^ [1^]. The right picture of fig.^ shows only solutions, which fulfill 

the dark matter condition. For bino-like LSP neutralinos the cosmologically allowed mass 
range for the lightest scalar Higgs 5*1 is 40 GeV < ms^ < lAOGeV, which follows from the 
cosmologically allowed mass range for the lighter selectron of 100 GeV < < 300 GeV^. 
Very light scalar Higgses 5*1 can be cosmologically allowed, if the LSP neutralino is a mixed 



singlino or a nearly pure singlino with mass m^o = so that Z-pole enhanced 

neutralino pair annihilation occurs. The corresponding plots in the MSSM [^, ^ look 



very similar with the exception, that there are no Higgses below the bound of 62.5 GeV 



v) Conclusions. We have investigated the particle spectra of the MSSM and NMSSM 
without and with the dark matter condition. In contrast to the MSSM, the squark 
masses in the NMSSM are nearly proportional to the gluino mass. We found, that pure 
LSP singlinos and most of the light scalar Higgs singlets are cosmologically excluded. 
With the dark matter condition the following cosmologically allowed mass ranges for 
the SUSY particles in the NMSSM can be given: m^.. < 300 GeV, m-^^ < 300 GeV, 
300 GeV < rriuj, < 1900 Gel^, 200GeV < mi^ < 1500 GeV, 350GeV < m-g < 2100 Gel^, 
ms^ < 140 GeV. The upper bounds for the mass of the SUSY particles in the NMSSM 
are comparable to the MSSM bounds. The larger Higgsino or singlino portion of the 
lightest LSP neutralinos in the NMSSM leads to improved lower bounds for the gluinos 
and squarks compared to the MSSM. 
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Figure Captions 

Fig 1 The mass of the LSP neutralino versus the lighter selectron mass. The upper 
pictures are for the MSSM, the lower pictures are for the NMSSM. In the right 
pictures the dark matter condition is imposed. 

Fig 2 The mass of the lighter u-squark and the lighter stop versus the gluino mass in the 
NMSSM. In the right picture the dark matter condition is imposed. 

Fig 3 The mass of the lightest Higgs Si versus the gluino mass in the NMSSM. In the 
right picture the dark matter condition is imposed. 
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Fig. 1 
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Fig. 2 
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Fig. 3 
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